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A unitary model of the final state Kn interaction amplitudes in the B — > Knn decays is con- 
structed. The weak decay penguin amplitudes, derived in QCD factorization, are supplemented 
by phenomenological contributions. The strange Kn scalar and vector form factors are used to 
calculate the Kn effective mass and helicity angle distributions, branching ratios, CP asymmetries 
and the phase difference between the B° and B° decay amplitudes to K*(892)n. The fit on the 
phenomenological parameters leads to a good agreement with the experimental data, particularly 
for the B -> K* (892)n decays. However, our predicted B ± -> K% (1430)7!;*, K% (1430) -> 
branching fraction is smaller than the results of the Belle and BaBar collaborations, obtained 
from isobar model analyses. A new parameterization of the 5-wave Kn effective mass distribu- 
tion, which can be used in future experimental Dalitz plot analyses, is proposed. 
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1. Introduction 

Studies of three-body charmless hadronic decays of B mesons are very useful not only in 
standard model tests and in searches for "new physics" effects but also in the determination of 
strong interaction amplitudes. Both weak and strong interactions can create structures seen on 
Dalitz plots. Analyses of these diagrams should be performed in a unitary approach which allows 
for a proper construction of B-decay amplitudes. Usually a suitable partial wave analysis of final 
state amplitudes should be done. Then an adequate determination of branching fractions and CP 
asymmetries for different quasi-two body decay reactions is possible. 

Construction of a fully unitary three-body strong interaction amplitude is, however, a difficult 
task. A first step towards this goal is to enforce two-body unitarity. Here we apply this concept 
to the Kn channel for the B — > Knn decays by using a unitary coupled channel model. We study 
the Kn amplitudes in the limited Kn effective mass range smaller than about 1.8 GeV. Our aim is 
to describe such physical quantities as differential effective mass and helicity angle distributions, 
integrated branching fractions and direct CP asymmetries. 

In experimental analyses of B-decays the isobar model is very frequently applied. Within 
that model quasi-two body branching fractions are determined. However, the decay amplitudes 
commonly used in the isobar model are not unitary neither in three-body decay channels nor in 
two-body subchannels. This lack of unitarity can create severe problems in the determination of 
branching fractions in the case of wide overlapping resonances. In B —> Knn decays one observes 
a wide S-wave resonance K$ (1430). Its width equals to about 270 MeV and the postulated K$ (800) 
state can have even larger width of 500 MeV. Thus an important source of model errors in extraction 
of the branching ratio for the decay B ± — ► (1430) 71 * is a possible wrong attribution of a part of 
Dalitz plot density to a background amplitude and to its interference with other amplitudes, mostly 
with the S-wave. Below, we shall briefly discuss that issue. 



2. Theoretical model 

The weak decay amplitudes of B ± , B° and B Q , which are derived in QCD factorization, are 
supplemented by phenomenological contributions to the penguin amplitudes. Strong interaction 
amplitudes are constrained by chiral symmetry, QCD and experimental data on meson-meson in- 
teractions. The matrix elements of the effective weak Hamiltonian involve the strange Kn scalar 
and vector form factors. The introduction of form factors, constrained by theory and other experi- 
ments than B decays, is an alternative to the use of the isobar model. The Kn 5-wave contribution 
to the B — > K~n + n~ amplitude reads: 
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Expressions for other contributions to the B decay amplitudes can be found in [jj]]. In the above 
equation Gf denotes the Fermi coupling constant, Mb, tuk and m n are the masses of the charged 
B mesons, kaons and pions, 71 (q 2 ) and fff K (q 2 ) are the B to %~ and the Kn scalar form 
factors. The symbols X u = V u bV* s ,X c = V c bV* s , are products of the Cabibbo-Kobayashi-Maskawa 
quark-mixing matrix elements V qq i; a"' e (S), j = 4,6,8, 10 are the coefficients of the effective Hamil- 
tonian H e ff, q 2 is the Kn effective mass squared and my, m s and are b, s and d quark masses. 
Finally the cl,c c 4 , c\ and c c 6 are the phenomenological complex parameters which are fitted to the 
experimental data. 

The scalar and vector Kn form factors are connected to scattering amplitudes in the S and 
P waves via unitarity relations. In the S wave we treat two coupled Kn and Ki]' amplitudes. 
Three coupled channels: Kn, K*n and Kp appear in the P wave. The corresponding scattering 
amplitudes are constrained by experimental data, especially by the LASS results obtained at SLAC. 
The Muskhelishvili-Omnes equations are used to calculate the Kn form factors. 

3. Results and discussion 

The Kn effective mass and helicity angle distributions, branching ratios, CP asymmetries and 
the phase difference between the B° and B° decay amplitudes to K* (892)71 are calculated using the 
minimization program which serves to fit the set of 319 data from the Belle and BaBar collabo- 
rations with the four phenomenological parameters c"'i. As input we use only the well measured 
branching fractions for the B — > K*(S92)n. We do not use the experimental branching fractions 
for the B — * ^(1430)71, which are not well determined due to the large width of the ^(1430) 
resonance. Some of our model predictions for the 5-wave part of the branching fraction for the 
B — ► Knn decays are given in Table 1 . One can notice that they are lower than the experimental 
data, being substantially smaller than the Belle results and closer to the BaBar numbers. If the 
parameters c^'g are all put equal to zero then the P-wave part of the branching fraction is underes- 
timated by a factor of 4 to 5 and the 5-wave part by a factor of 2. The fit on the model parameters 
leads to a good agreement with the experimental data, particularly for the kaon-pion effective mass 
and helicity angle distributions. Some results are presented in Fig.l. 

Based on a good description of data we propose to choose in the future experimental analyses 
the following parameterization of the 5-wave B — > Knn amplitude: 

=fo*(*4x)(po/n% x + ci), (3-D 

where is the scalar Kn form factor while cq and c\ are complex numbers to be fitted from the 
data. Numerical values of the complex scalar form factor can be provided on request. 
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Table 1: Branching ratios averaged over charge conjugate reactions B — > Kkk in units of 10~ 6 . {K + n~)p 
and [K + n~)s denote the (K + n~) pair in f-wave and S-wave, respectively. In the first two lines, the values 
of the model, calculated by the integration over the given rriKn range, are compared to the corresponding 
Belle [fi ||] and BaBar [0, [|] results written in the third and fourth column. In the last two lines, the Belle 
branching fractions [^, calculated with a ^(1430) Breit-Wigner amplitude, and the BaBar branching 
fractions [|[ |J], calculated in their parametrization of the Kn 5-wave, are compared to our model predictions. 
The model errors are the phenomenological parameter uncertainties found in the minimization procedure. 



decay mode 


m Kn range (GeV) 


Belle 


BaBar 


model 


B+ -» {K+K-) P 1T 


(0.82,0.97) 


5.35 ±0.59 


5.98 ±0.75 


5.73±0.14 


B° -» (K°n + ) P %- 


(0.82,0.97) 


4.65 ±0.77 


6.47 ±0.75 


5.42±0.16 


B+ -» (K+k-) s K- 


(0.64,1.76) 


27.0 ±2.5 


22.5 ±4.6 


16.5±0.8 


B° -► (K Q n+) s %- 


(0.64,1.76) 


26.0 ±3.4 


17.3±4.6 


15.8±0.7 



200 



O 



100 



CO 

+-> 

d 

> 
CD 



50 







" 1_r T 

a) 



M 1 1 1 i 1 1 1 i 1 1 1 i 1 

B* "> K* TT T 7T ± 




0.6 0.8 1.0 1.2 1.4 

( GeV ) 



1.6 




Figure 1: a) The K 71^ effective mass distributions in the B — > K^if^it decays. Data points are from 
Ref. [Q], The dashed line represents the 5- wave contribution of our model, the dotted line that of the P- 
wave and the histogram corresponds to the coherent sum of the S- and P-wave contributions, b) Helicity 
angle distribution for B — > K ± n :f n ± decays calculated from the averaged double differential distribution 
integrated over m K ± n ^ mass from 1.0 to 1.76 GeV. Data points are from Ref. [H], The dashed line represents 
the S-wave contribution of our model, the dotted line that of the P-wave and the dot-dashed line that of the 
interference term. The histogram corresponds to the sum of these three contributions. 
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